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Abstract 

Scattering from a compound barrier, one composed of a number of distinct 
non-overlapping sub-barriers, has a number of interesting and subtle mathe- 
matical features. If one is scattering classical particles, where the wave aspects 
of the particle can be ignored, the transmission probability of the compound 
barrier is simply given by the product of the transmission probabilities of the 
individual sub-barriers. In contrast if one is scattering waves (whether we are 
dealing with either purely classical waves or quantum Schrodinger wavefunc- 
tions) each sub-barrier contributes phase information (as well as a transmission 
probability), and these phases can lead to either constructive or destructive 
interference, with the transmission probability oscillating between nontrivial 
upper and lower bounds. In this article we shall study these upper and lower 
bounds in some detail, and also derive bounds on the closely related process 
of quantum excitation (particle production) via parametric resonance. 

Keywords: Compound barrier, transfer matrix, transmission amplitude, 
transmission probability, reflection amplitude, reflection probability, Bogoli- 
ubov coefficients, quantum particle production, parametric resonance. 
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1 Background 



Consider a general one- dimensional scattering problem. One might be interested in 
classical waves (acoustic waves in a pipe, surface waves in a channel, electromagnetic 
waves in a waveguide), or quantum waves (the Schrodinger equation). One formalism 
that is very well adapted to addressing this mathematical situation is that of "transfer 
matrices" where one relates the waves on the left of a barrier to the waves on the right 
of the barrier via a 2 x 2 complex transfer matrix. These transfer matrix techniques 
are discussed, at varying levels of detail, in the textbooks by Merzbacher [TJ, Mathews 
and Venkatesan [2], and Singh |3j, and in the pedagogical articles by Sanchez-Soto et 
al [I] and the present authors |5J. Other more technical research articles using this 
formalism include those of Peres [6], Kowalski and Fry [7], Korasani and Adibi [8], 
and Barriuso et al [9j [TO] • These techniques have also served as backdrop to some 
previous results reported by the second author in [IT] , and related work by both 
current authors reported in^2l[T3l[Tl[T5l[T6l[T71[T8]. 

We start by noting that transfer matrices are of the form: 

«| 2 -|/3| 2 = 1, (1) 

where a and (3 are known as "Bogoliubov coefficients" . In terms of the perhaps more 
common transmission and reflection amplitudes one has 

t| 2 + |r| 2 = l. (2) 

The transmission and reflection probabilities are simply 

T=\t\ 2 ; R=\r\ 2 ] T + R=l. (3) 

Because of these normalization properties we can always find real numbers {9, <p a , <pp} 
to write 

a = cosh(0) exp(i(p a ); (3 = sinh(0) exp(i<pp); (4) 

and 

t = sech(#) exp(— i<p a ); r = tanh(6 l ) exp(— i[<p a — ipp]). (5) 

Now consider a situation where one has a number of non-overlapping barriers. We 
can assign a distinct transfer matrix M, to each barrier. How do these transfer 
matrices combine to give the "total" transfer matrix for the compound barrier? 



M 



a (3 
(3* a* 



M 



l/t r/t 
r*/t* l/t* 
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Two key standard results are [U EJ [3J HJ [5] : 



If one moves any one of the sub-barriers a distance Oj, the effect of this trans- 
lation is to modify 



Mi -> M 



'rikai 









Mi 



-ikcii 







i-ikai 



on f3 ie +2lka > 
f3*e- 2ika > a* 



(6) 

That is, the Bogoliubov coefficient (and so the transmission amplitude tj) 
is invariant under a shift in the physical location of the barrier, while the 
Bogoliubov coefficient /3j (and so the reflection amplitude ri) picks up a shift- 
dependent phase e +lkai . 

• A compound transfer matrix (for n distinct non-overlapping localized barriers) 
is of the form 

M 12 ... n = M X M 2 ...M n . (7) 
The order in which these matrices are multiplied together is important. 

The various phases e %LPa ^ , e i</3/3 » , and e +tkai , have a very strong influence on the overall 
transmission and reflection properties, and can lead to either constructive interfer- 
ence or destructive interference. (In extreme cases one may encounter transmission 
resonances [reflection-less potentials] where the transmission probability is unity, 
T — > 1, even for rather nontrivial potentials.) 

We shall be interested in extracting as much information as possible regarding 
the compound transmission probability T 12 ... n without having access to knowledge 
of the individual phases e ¥ "' , e i</3/3 *, and e +tka \ That is, given only the individual 
transmission probabilities Tj, what can one say about the total transmission prob- 
ability Ti2...n? Naturally, information regarding transmission probabilities will also 
translate into information regarding reflection probabilities. 

We have so far phrased things in terms of a scattering problem, but if one moves 
from the space domain into the time domain, the same analysis will give information 
regarding quantum particle production due to parametric resonance. Each "barrier" 
is now viewed as a time- dependent parametric change in oscillation frequency, and 
each of these "parametric excitation episodes" would by itself excite a number iVj of 
quanta, where 

Ni = m 2 . (8) 

But, because of the phases e tipa ^, e^, and e +lkai , distinct parametric excitation 
episodes can interfere either constructively or destructively, and the total number of 
quanta excited, iV 12 ... n , will be some complicated function of the individual Ni and 
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the phases. (In extreme cases one may encounter total destructive interference, where 
the nett (overall) particle production is zero, iV — > 0, even for a rather nontrivial 
parametric excitation.) 

We shall also be interested in extracting as much information as possible regard- 
ing the nett particle production Ni 2 ... n without having access to knowledge of the 
individual phases e i</3 % e 1 ^*, and e +tkai . That is, given only the individual particle 
production numbers iVj, what can one say about the total particle production N\ 2 _ n l 

In counterpoint to these phase-dependent questions, if one is working in the limit 
where the objects being scattered can be treated as classical particles, so their wave 
properties can be neglected, then the individual transmission probabilities Tj are 
combined according to the standard rules of ordinary probability theory: 

n 

-^particles I I T%- (9) 
i=l 

In the discussion below we will investigate these questions in considerable detail, 
developing very general upper and lower bounds on the nett Bogoliubov coefficients, 
nett transmission and reflection probabilities, and nett number of excited quanta in 
terms of the properties of the individual sub-barriers and/or parametric resonance 
episodes. 

We have also investigated the possibility of utilizing the Abeles matrix formula- 
tion of the scattering problem. (See for instance [T91 I20~t [2TJ [22].) While the transfer 
matrix approach and the Abeles matrix approach are fundamentally just different 
mathematical ways of addressing the same physical problem, and so must ultimately 
give completely equivalent answers, the the details of the matrix manipulations in- 
volved are different. Specifically, the process of extracting a transmission amplitude 
from an Abeles matrix is somewhat messier than that for extracting a transmission 
amplitude from a transfer matrix — for this reason we focus on the transfer matrix 
approach throughout the current article. 

2 Two-barrier systems 

The compound transfer matrix (for two non-overlapping localized barriers) is of the 
form 

ai« 2 + Pi(3 2 ai/3 2 + (3ia 2 
alPZ + Pta 2 ala* 2 + (3l(3 2 



M 12 = M 1 M 2 

That is: 



(10) 

«i2 = oi\a 2 + f3i/3 2 ; f3 12 = ai/3 2 + (3ia* 2 . (11) 
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2.1 Bounding the Bogoliubov coefficients 

This leads to the immediate upper bounds 

|a 12 | < | ai ||a 2 | + |/?i 2 | < |«i||/3 2 | + |/3i||« 2 |. (12) 

Now using the normalization conditions, and writing = cosh#j and — sinh^j, 
we have 

|a 12 | < cosh #x cosh # 2 + sinh^i sinh# 2 = cosh(6*i + 6> 2 ), (13) 

and 

\/3 l2 \ < cosh 6»i sinh 6» 2 + sinh 6»i cosh 6> 2 = sinh(6>i + 8 2 ), (14) 
which can be rewritten as 

1 0:12 1 < coshjcosh -1 I q;x I + cosh -1 | ct 2 1 } ; (15) 

and 

\(3 12 \ < sinhjsinh- 1 |^| + sinh" 1 \(3 2 \}. (16) 
In the other direction we have the immediate lower bound 

K2I > KIM - I /5i 1 1 /5 2 1 , (17) 

which we can rewrite as 

I cux2 1 > coshjcosh -1 \cti\ — cosh -1 \a 2 \}. (18) 

Because cosh is an even function of its argument, there is no need to worry about 
the relative sizes of the two cosh -1 \ai\ terms. In contrast, we have to be particularly 
careful when deriving a lower bound for /3 12 . (This will be a recurring theme in the 
article, so we shall be particularly explicit in deriving this result.) Note that 

p 12 = ai p2 + Pict* 2 = |cki||/3 2 | exp (i[<p ai + (ffo]) + \Pi\\a 2 \ exp(i[^ 1 - ip a2 ]) . (19) 
But then 

IA2I = \ai\\/3 2 \ + |/3i||o!2|exp(i[^ / 3 1 - tp^ - <p ai - <p a2 ]) . (20) 

But I o; 1 1 1 /5 2 1 and |/3i||ai2| are both real and positive - - so the quantity above is 
minimized when destructive interference is maximum. It is now easy to see that this 
occurs when the phases are such that 

<ppi - - <Pai - = ( 2n + iK (2i) 
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in which case 



exp(t 



Consequently 



1 > KHAd - \Pi\\a 2 \ 



(22) 
(23) 



Note that the outermost set of absolute value signs is particularly critical, guaran- 
teeing that the RHS is non-negative. This can now be rewritten as 



|/5i2 1 > sinh sinh" 1 - sinh" 1 \(3 2 \ 



(24) 



Again, the outermost set of absolute value signs is particularly critical, guaranteeing 
that the RHS is non-negative. 

2.2 Bounding transmission and reflection probabilities 

Now working in terms of transmission and reflection amplitudes and probabilities we 
have 

(25) 
(26) 



M 



a (3 




- 1/t 


r/t 


(3* a* 




r*/t* 


1/t* 



and 



\t\ 



\a\ 



sech 2 9; R 



\a\ 



tanh 2 6. 



The upper bounds on the Bogoliubov coefficients lead to a lower bound on T and an 
upper bound on R as follows: 



T 12 > sech 2 |sech" 1 v / T\ + sech _1 v / ^} ; 
R 12 < tanh 2 jtanh" 1 ^/r[ + tanh" 1 • 



(27) 
(28) 



Similarly, the lower bounds on the Bogoliubov coefficients lead to an upper bound 
on T and a lower bound on R as follows: 

(29) 
(30) 

Note that, because we are squaring the sech and tanh, we do not need to worry about 
the relative magnitudes of the Tj and the Ri in the two formulae above. 



T12 < sech 2 | sech — sech 1 a/7^| ; 

Rn > tanh 2 jtanh -1 \/R~i — tanh -1 \/i?2 j ' 
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By manipulating the hyperbolic and inverse hyperbolic functions, (see some key 
hyperbolic identities in appendix EJ, these bounds can be brought into the rational 
algebraic form 



Ti 2 > 



T{T 2 



{1 + yJT=T iy /T=%} 



2 ' 



#12 < 



Similarly 



Ti 2 < 



T X T 2 



{1- V^TiV^T^} 



2 ' 



#12 > 



y/R~l + V#2 

1 + v^v 7 ^ 



^R~i — \fll2 



(31) 



(32) 



As a useful internal consistency check on the formalism note that for Ri £ [0, 1] we 
have 

a ' #1 4" v #2 



and that for Tj e [0, 1] we have 



1 + VR1VR2 



TiT 2 



< 1, 



{i-yr^yr^K 



< 1. 



(34) 



Note that to get an exact transmission resonance Tj_2 = 1 it is necessary (but certainly 
not sufficient) that 

™ s = 1. (35) 

{1- v/i^vr^} 2 

But this implies Ti = T 2 . That is, if one ever wishes to obtain an exact transmission 
resonance from two disjoint non-overlapping barriers, then the two individual barriers 
must have equal transmission probabilities. 



2.3 Bounding particle production 

If in contrast we work in a particle production scenario, (via episodic parametric 
resonance), where N = \(3\ 2 we can similarly extract upper and lower bounds 



iVi2 < sinh 2 {sinh" 1 \fN~i + sinh" 1 v 7 ^} 
N n > sinh 2 jsinh" 1 y/Wx - sinh" 1 1/^2} 



(36) 
(37) 
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Again, because we are squaring the sinh, we do not need to worry about the relative 
magnitudes of the Ni in the formula above. These bounds can be converted (see 
appendix [A} to the algebraic statements 

N 12 < {y/Ni(N 2 + 1) + ^/N^Nt + 1)}' ; (38) 

N l2 > {y/N^ + 1) - y/N^h + 1)}' • (39) 

Note that if one wishes to get complete destructive interference, then one must 
have ^■ s /Ni(N 2 + l) - \/ N 2 (N 1 + 1) j = 0, implying JVi = N 2 . That is, if one ever 
wishes to obtain an exact cancellation of particle production from two disjoint non- 
overlapping parametric resonance episodes, then the two individual episodes must 
individually have equal quantities of particle production. 



2.4 Summary 

Having worked through the two barrier case in some detail, we are now in a position 
to consider the more general case. Note that all the bounds above have been care- 
fully checked against the extant literature, (see for example [Tj, [3J, SI El El El El El 
HDl HD H21 U21 El 151 US HO Hg), and are compatible with all known results. For 
the multiple-barrier situation, many of the calculations are immediate and straight- 
forward generalizations of the above — we shall soon see that the tricky one is the 
lower bound on the (absolute magnitude of the) Bogoliubov coefficients |ai2...n| and 
|/3i2...n|, (which ultimately leads to the upper bound on T\ 2 .„ n and the lower bound 
on N 12 ___ n ). 



3 General upper bounds on \a\ and \/3\ 

The two-barrier upper bounds on the Bogoliubov coefficients \a\ and \{3\ immediately 
generalize (by straightforward iteration) to the composition of n transfer matrices 

I«i2...n| < cosh jy^cosh" 1 \aj\\ ; \Pi2...n\ < sinh | sinh" 1 |^|| . (40) 

This immediately leads to a lower bound on transmission probability and an upper 
bound on reflection probability. 

T 12 ... n > sech 2 | ^ sech_1 V^j ! Ri2...n < tanh 2 j tanh" 1 ^/r\ . (41) 
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If we work in a particle production scenario where N = \(3\ 2 we can similarly write 

N 12 ... n < sinh 2 | sinh_1 J . (42) 

In the general case these appear to be the most tractable form of the bounds. (We 
have also verified the correctness of these results by explicitly iterating the two-barrier 
results to investigate three-barrier, four-barrier, and certain particularly tractable n- 
barrier systems, comparing with the general result presented above.) 

4 General lower bounds on \a\ and \/3\ 

When it comes to bounding the Bogoliubov coefficients from below we have already 
seen 

I«i2| > |«i||« 2 | - lAH&l; IA2I > \<xi\\p2\ ~ \/3i\\®2\ ; (43) 
which can be rewritten as 

I 1 > cosh{cosh _1 \cti\ — cosh" 1 \a 2 \}] \/3 12 \ > sinh sinh -1 — sinh" 1 \/3 2 \ ■ 

(44) 

However generalizing these lower bound inequalities to n steps is much more difficult. 
We shall build up our lower bound in stages, first providing a recursive version of 
the bound before constructing an explicit solution to the recursion relation. 

4.1 Implicit iterative lower bounds on \a\ and \(3\ 
Define a collection of n parameters: 

9i = cosh -1 I a; I, (45) 

and the sums (m G {1, 2, 3, . . . n}) 

m 

S m = Y, 9 *- (46) 
i=i 

Now take 

B 1 = 9 1 , (47) 
and, for m G {1, 2, 3, ... n — 1}, iteratively define the quantities B m+ i by 

B m +i = (Om+i ~ S m ) H{O m+ i — S rn ) + (B m — Orn+i) H{B rn — 6 m +i), (48) 
where if (•) is the Heaviside step function. 
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Theorem: By iterating the 2-step bounds one has 

| > coshS n ; \Pi2.. , n \ > sinh_B n ; (49) 

and consequently 

T 12 ... n < sech 2 5 n ; R 12 ,„ n > tanh 2 B n ; N 12 ___ n > sinh 2 B n . (50) 

Proof by induction: When we iterate the definition for B n the first two times we 
obtain 

51 = 61, B 1 = 6 1 ; (51) 

5 2 = 9 1 + 9 2 ; B 2 = \B l -9 2 \. (52) 

Thus by our previous results, the claimed theorem is certainly true for n = 2. Now 
apply mathematical induction: Assume that at each stage the interval [B m , S m ] 
characterizes the highest possible and lowest possible values of 6i 2 ... m . Applying the 
2-step bound to the pair 0\2...m an d 6 m +i leads trivially to 6> 12 ...( m+ i) being bounded 
from above by 

S m +i = S m + m +i, (53) 
and less trivially to being bounded from below by 

B m+ i = (6 m +i — S m ) H(9 m+ i — S m ) + (B m — 9 rn+ i) H{B rn — 6 m +i), (54) 
as asserted above. This completes the inductive step. That is: 

012...(m+l) ^ [B m +i, S m+ i\. (55) 

□ 

However these bounds are unfortunately defined in a relatively messy and implicit 
iterative manner — can this be usefully simplified? Can we make the bounds explicit? 

4.2 Symmetry properties for the lower bound 

When we iterate the general definition of B n , previously used to obtain B 2 , one more 
time we see 

s 3 = e 1 + e 2 + e 3 - B 3 = max{e 3 -(9 1 + e 2 ),\9 1 -e 2 \-e 3 ,o}. 

Is there any further way of simplifying this? Rewrite B 3 as 

B 3 = max{9 1 - 6 2 - 9 3 , 9 2 - 9 3 - 6 U 9 3 - 9 X - 9 2 , 0}. (56) 

Note that this form of B 3 is manifestly symmetric under arbitrary permutations of 
the labels 123. One suspects that there is a good reason for this. In fact there is. 
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Theorem: Vn B n {9i) is a totally symmetric function of the n parameters 0j. 

Proof: Note that the individual transfer matrices Mj and Mj trivially have the 
same values of Si and B{, and in fact have the same values of 0j. Note further that 
for any two transfer matrices MiM 2 and M 2 Mi have identical values of S 2 and B2. 
Finally for any two transfer matrices the products M\M 2 and M\M 2 have identical 
values of S 2 and B 2 . (These assertions all follow from the simple results obtained 
above for compounding two transfer matrices.) 
That is: 

S(Mi) = S(MT); £(M) = B(M?); (57) 
S{M 1 M i ) = S(M 2 M 1 ); S(MiM 2 ) = B(M 2 M 1 ); (58) 
,S(MiM 2 T ) = 5(MiM 2 ); B(M 1 M 2 ) = B{M 1 M 2 ). (59) 
But then by cyclic permutation 

B(M 1 M 2 M 3 ) = B(M 3 M 1 M 2 ) = B(M 2 M 3 M 1 ), (60) 

and furthermore 

B(M 2 M 1 M 3 ) = B((M^M^) T M 3 ) 
= B((M± M 2 )M 3 ) 
= B(Mj (M 2 M 3 )) 
= B(M 1 (MjM 3 )) 
= B((M 2 M 3 )Mi) 
= B(M 2 T (M 3 Mi)) 
= B(M 2 M 3 M 1 ) 

= B(M 1 M 2 M 3 ). (61) 

That is 

B(M 2 M 1 M 3 ) = B(M 1 M 2 M 3 ). (62) 

Combining these results implies that B(MiM 2 M 3 ) is a symmetric function of the 
three transfer matrices Mj, and hence a symmetric function of the three parameters 
0i, 02, 0,3- But this argument now generalizes — For any B{M\M 2 . . . M n ) you can 
use this argument to show 

B(M 1 M 2 M 3 ...M n ) = B(M 2 M X M 3 . . . M n ); (63) 
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and 

B(M 1 M 2 M 3 ...M n ) = B(M 2 M 3 . . . M n Mi), (64) 

which implies complete symmetry in all its n arguments 0j. □ 
Based on this observation we can now assert a bolder theorem that has the effect of 
yielding an explicit (non-iterative) formula for B n . 

4.3 Explicit non-iterative formula for the lower bound 

Theorem: 

Vm : B m = max {29 { - S m , 0}. (65) 

ie{l,2,...m} 

Proof by induction: We have already seen that the iterative definition of B n can 
be written as 

B m+ i = max{# m+ i — S rn , B rn — 9 m+1 , 0}, (66) 
which we can also rewrite as 

B m+ i = max{26' m+ i — S m +i, B rn — 9 m +i, 0}. (67) 

Now apply induction. The assertion of the theorem is certainly true for m = 1 and 
m = 2, and has even been explicitly verified for m — 3. Now assume it holds up to 
some m, then 

B m +i = max{2^ m+1 — S m+ i, B rn — 9 m+ i, 0} 

= max <^ 29 m+1 - S m+1 , max {26 i - S m , 0} - 9 m+1 , 

[ ie{l,2,...m} 



= max <^ 29 m+1 - S m+1 , max {26 { - S m+1 , 0}, 

I. i6{l,2,...m} 

= max {2^-^+1,0}. 

ie{l,2,...m,(m+l)} 



This proves the inductive step. Consequently 



and in particular 



Vm : B m = max {26 i — S m , 0}, (68) 

i6{l,2,...m} 



B n = max {29i-S n ,0}. (69) 

ie{l,2,...n} 
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(For completeness, note that we have explicitly checked the equivalence of these 
iterative and non-iterative results for B n by using symbolic manipulation packages 
for multiple examples for the cases n = 5 and n = 10.) □ 
To simplify the formalism even further, let us now define 

# P eak = max $ h (70) 

ie{l,2,...n} 

and 

n 

^offpeak = 0i = 2_j ® l ~ ^P cak = ^ n _ ^Peak- (71) 

(We wish to use the subscript "peak" for the maximum of the individual 0j's; the 
subscripts "max" and "min" will be reserved for bounds on the n-fold composite 
transfer matrix.) Then we can write 

B n = max{26 l p 0ak - S n , 0}, (72) 

or 

B„ = max{6 l p ea k — o ffpeak, 0}. (73) 

Note that the max function is still needed to guarantee that the B n are non-negative. 
With this explicit formula for B n in hand, we have 

I«i2...n| > cosh [max{26 l p eak - S n , 0}] ; \/3 12 ... n \ > sinh [max{26» peak - S n , 0}] . 

(74) 

4.4 Geometrical interpretation 

These lower bounds on the Bogoliubov coefficients can also be given a clean and 
intuitive geometrical interpretation. From the work of Barriuso, Monzon, Sanchez- 
Soto, and Carinena [U EJ [10] it is known that the composition of two barriers can be 
understood as the composition of two hyperbolic vectors (translations in the hyper- 
bolic plane, "turns" in Hamilton's original notation). The length of these hyperbolic 
vectors is related to the Bogoliubov, transmission, and reflection coefficients by 

£i = cosh -1 \ai\ = sinh -1 = sech _1 -\/7^ = tanh -1 \^Ri. (75) 

The two-barrier upper bounds ( Tl3l) and (fl4l) can then be easily derived from geo- 
metrical considerations as hyperbolic space applications of the triangle inequalities. 
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The lower bounds ( IT7l) and ( 1241) are a little more subtle, but for the two-barrier case 
both upper and lower bounds can be summarized as 



\k-l2\<ln<k+k- (76) 
For the n-barrier case the upper bound is straightforward 



-12. ..n 

i=l 



The lower bound is a little trickier. Separate out the lengths li of the hyperbolic 
vectors into the "largest" and the "rest" (corresponding to what we previously called 
"peak" and "off peak"). Then geometrically 

> max <^ £i ar g CS t - ^Ji, > . (78) 

L rest J 



-12.. .n 



It is only when the single largest step exceeds the maximum possible size of all 
the other remaining steps put together that one obtains a non-trivial lower bound. 
We shall now apply this formalism to bounding the the transmission and reflection 
probabilities, and to bounding the amount of particle production. 

5 Application to T, i?, and N 

First, we note that 

Tl2... n < Sech 2 {max{6> p eak - #off P eak, 0}} , (79) 

or equivalently 

T 12 ... n < sech 2 {max{2# pcak - S n , 0}} . (80) 

That is 

Tn...n < sech 2 |max|2 sech -1 ^/T pcak - sech -1 A/T min , o| | , (81) 
where we have defined T min by 

Tvi...n > T min = sech 2 {S n } . (82) 

Second, we note that 

R\2...n > tanh 2 {max{26' peak - S n , 0}} , (83) 



15 



that is 

R\2...n > tanh 2 | max |2 tanh" 1 ^ 'R pcak - tanh -1 a/ i? max , o| j , (84) 
where we have defined -R max by 

Ri2... n < ^max = tanh 2 {S n } . (85) 

Finally, 

N 12 ... n > sinh 2 {max{2# pcak - S n , 0}} , (86) 

that is 

N 12 ... n > sinh 2 {max |2 sinh" 1 y/N peeik - sinh' 1 v^max, o} } , (87) 
where we have defined iV max by 

N 12 ...n < ^max = sinh 2 {S n } . (88) 

If we want to be a little more explicit, then we could write these bounds as 

T 12 ... n < sech 2 jmax 1 2 secli -1 y/T^ - sech ^ V^i, oj J ; (89) 

R\2...n > tanh 2 J max J 2 tanh -1 v^R^ - ^ tanh -1 1 1 ; (90) 



i=i 



A^i2...n > sinh 2 jmax 1 2 sinh" 1 - ^ sinh" 1 v 7 ^, j j . (91) 

This appears to be the simplest form achievable for these three bounds. Remember 
that these three bounds, coming from lower bounds on the Bogoliubov coefficients 
are complemented by three other simpler bounds 

T 12 ... n > sech 2 | sech_1 V^j 5 Ri2...n < tanh 2 | tanh -1 v^j ; (92) 

N 12 ... n < sinh 2 | s^ 11 ' 1 V^j • (93) 
coming from upper bounds on the Bogoliubov coefficients. 
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6 Perfect transmission resonances 



Transmission resonances (perfect transmission resonances) occur when T — > 1 at 
certain energies or barrier spacings when the phases work out just right to cancel the 
reflection. But the occurrence of these resonances is still constrained by our general 
bound 

T\2...n < sech < max |2sech- 1 v ^^-^sech~ 1 v / ^,o||. (94) 
So a perfect transmission resonance can only occur if 

n 

2sech- 1 v ^^-^sech- 1 v^< 0, (95) 
i=i 

(this is a necessary condition, not a sufficient condition). This is equivalent to re- 
quiring 



T P eak < sech 2 I ^secli- 1 ^! = sech 2 jisech- 1 



(96) 



Now use the "half angle formula" 



sech(x/2) = J— -= = J 2seGh f , (97) 

J coshrr + 1 V 1 + sechx' V ; 



to obtain 



Again, this is a necessary condition (not a sufficient condition) in order for a perfect 
transmission resonance to be possible. 

In terms of a particle production scenario, a transmission resonance translates 
to "complete destructive interference" between two or more particle creation events 
(so that the nett particle production is zero). To be absolutely certain of avoiding 
"complete destructive interference" one must have 



> sinh 2 < - sinh 1 JW; \ . 



iVpeak > sinh 2 { - } ^ sinh" 1 ^N t } . (99 ) 

But we know 

N < N max ee sinh 2 S n = sinh 2 <j sinh" 1 \ , (100) 



i=i 
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so we can rewrite this as 



Npeak > ~ • (101) 

This is a sufficient condition for iV min > 0, a sufficient condition for there to be at 
least some overall particle production. 



7 Discussion 

We have seen how a rather nontrivial application of the transfer matrix formal- 
ism, (see for example [TJ |2l [3J SJ [5] for basic background information), allows us 
to place rather nontrivial bounds on the physics of compound scattering systems 
- even in the absence of detailed information concerning their internal structure. 
The compound scattering systems we are interested in are one- dimensional systems 
built up from a number of disjoint non-overlapping barriers; for such systems we 
have explicitly shown that the nett scattering properties (unsurprisingly) depend 
not only on the transmission probability, Tj, of each individual barrier, but also 
on a number of phases associated with these individual barriers and the separa- 
tion between them. Surprisingly, even in the absence of quantitative information 
regarding the value of these phases it is nevertheless possible to place rigorous and 
nontrivial upper and lower bounds on the scattering properties of the compound 
system. The resulting bounds also apply (with suitable modification of the lan- 
guage) to bounds on the number of quanta excited via parametric resonance. These 
bounds have all been carefully checked against the extant literature, (see for exam- 
ple PElEllllElElEllHlEll^IIIlITSlIIalllll^ QS]), and are compatible 
with all known results. The bounds are explicit, compact, and (despite their relative 
simplicity) appear to be entirely novel. 
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A Hyperbolic identities 

Here are several useful and quite elementary hyperbolic identities (particularly rele- 
vant to sections [2721 and [273]) : 

sinh (sinh -1 A + sinh -1 i?) = sinh (sinh -1 A) cosh (sinh -1 £?) 

+ cosh (sinh -1 A) sinh (sinh -1 £>) 
= A Vl + B 2 + Vl + A 2 B. (102) 

cosh (sinh -1 A + sinh -1 B} = cosh (sinh -1 A) cosh (sinh -1 i?) 

+ sinh (sinh -1 A) sinh (sinh -1 £>) 

= VTTa 2 VTTb 2 + A B. (103) 

cosh (cosh -1 A + cosh -1 = cosh (cosh -1 A} cosh (cosh -1 B^ 

+ sinh (cosh -1 Aj sinh (cosh -1 B) 
= AB + VA 2 - 1 VB 2 - 1. (104) 

tanh (tanh -1 A) + tanh (tanh -1 Z?) 



tanh (tanh -1 A + tanh 1 fi) 



1 + AB 

This now implies the more subtle result 

1 



1 + tanh (tanh 1 A) tanh (tanh 1 B) 

A + B , . 

(105) 



sech (sech 1 A + sech 1 B) 



cosh (cosh -1 (l/A) + cosh _1 (l/5)) 
1 



A^B- 1 + VA' 2 - WB- 2 - 1 

AB , . 
/ 106 
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